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We explore physics of unstable particles when mother particle mass is around the sum of its 
daughter particle masses. In this case, the conventional wave function renormalization factor is 
ill-defined. We propose a simple resolution of the threshold singularity problem which still allows 
the use of narrow width approximation by dehning branching ratio in terms of spectral density. The 
resonance peak and shape is different for different decay channels and no single decay width can be 
assigned to the unstable particles. Non-exponential decay happens in all time scales. 


Introduction The narrow width approximation (NWA) 
has played an important role in studying unstable par¬ 
ticles. Unstable particle states can not be asymptotic 
states of a scattering amplitude in order to keep unitar- 
ity and causality[1]. Nevertheless, NWA allows similar 
treatments of unstable particles by factorizing full scat¬ 
tering cross sections of stable states into production and 
decay parts. In most practical situations, heavy off-shell 
calculations are immensely simplified with NWA. 

When NWA is used for Standard Model calculations 
with realistic parameters, it is enough to take conven¬ 
tional wave function renormalization factor, Z for un¬ 
stable particles, whose inverse is defined by a real part 
of G~^ differentiated by momentum square at physical 
mass square. Conventional choice of the physical mass is 
a zero of Re(G“^). They work when all the dressed prop¬ 
agators [2] are well approximated by Breit-Wigner(BW) 
distribution [3]. 

However, this Z is ill-defined in some examples be¬ 
yond the Standard Model. In many cases, a self en¬ 
ergy included in a dressed propagator is proportional 

to = y^(l - ) where ma,b are 

masses of particles propagating in the loop. This is be¬ 
cause the phase space volume of decay is proportional 
to P and self energy and decay rate are closely related 
by the optical theorem. The classihcation of interactions 
providing the self energy with the same property is done 
later in this Letter. For simplicity, we discuss the prob¬ 
lem with a scalar theory in the text. 

For the self energy proportional to p, the Z~^ contains 
a term proportional to 1//3 which diverges as the physical 
mass approaches the threshold mass, ma-\-mb from below. 
Taken faithfully, Z —> 0 means every production and 
decay of the unstable particle vanish and the particle 
becomes isolated from the theory no matter how strong 
the interaction is, which is nonsense. 

Solving the ill-defined Z problem has been attempted 
mostly by using complex pole scheme[4, 5] which relates 
complex pole(s) on the second Riemann sheet to physi¬ 
cal quantities: its real part to physical mass, imaginary 
part to decay rate, and residue to Z. After the complex 
pole was conjectured to have physical meaning [6], its 


gauge independence was shown in Z boson in the Stan¬ 
dard Model [7] and scalars [8], and the scheme was em¬ 
ployed to Higgs physics[9]. However, using complex pole 
can be traced back to residue theorem for contour integral 
over lower half plane of the second Riemann sheet where, 
below a threshold, the analytically continued propagator 
G 2 dehned in the second Riemann sheet deviates from 
the correct propagator G which should have been used 
in exact calculation. 


To understand the problem, it is important to know 
what really happens to the dressed propagator as the 
physical mass approaches a threshold mass. If the phys¬ 
ical mass is near the threshold, the kinetic term can be 
represented by — m^ oc /3^, while the self energy term is 
proportional to p around the peak. Thus the self energy 
is dominant near the peak and it changes the shape of the 
propagater to be totally different from BW distribution. 
We propose generalized narrow width approximation by 
defining branching ratio in terms of spectral density. 


As the propagator G changes, p{p^), the spectral 

density of Kallen-Lehmann representation [10, 11] also 

changes. Unlike BW distribution which gives exponential 

decay with a rate of imaginary part of a complex pole, the 
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survival probability, P{t) = fp° p{S) does 

not exponential decay. Deviation from exponential de¬ 
cay for very short or long time in quantum held theory is 
well known [12, 13]. We show non-exponential decay pat¬ 
tern in middle range of time when most decay happens 
if they are at the threshold. 


Factorization Consider a full scattering cross section 
constructed with all external states by stable particles. 
For simplicity, assume that one Feynmann diagram (Fig. 
1) dominantly determines the process which contains 
an unstable particle state, p that ends up being a h- 
nal state A. After inserting the identity, / dSS{S — 
p's) f ~ P>P^{P%) iiito the full scattering cross 

section, we obtain 


cr( initial —> 1, 2, • • • , n. A) 

= / dS'(T(initial—> mi, •• • ,mn,'/S) p\{S) (1) 
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FIG. 1. Feynman diagram of initial —>■ 1, 2, • • • ,n,X. 

where A'max (•S'min) is the biggest (smallest) possible in¬ 
variant mass squared of A and 

Pa(5)^|(0|</)(0)|A;5)P (2) 

which is a spectral density to a specific channel A. If A 
contains more than one particle, V'S'max is — 
where Etot is total initial energy. When n > 1, we can 
take SjYiiYi — 0- 

If p\{S) is delta function like with center m\ and 
width mAhA, and the production part, CT(initial —> 
mi,-- - ,mn,V~S) does not change rapidly within the 
width, then we can approximate eq.(l) by ct( initial —> 
mi,-- - ,mn,rn\) dSp\(S). For example, with mo- 
mentum independent interaction, ^(initial —>■ mi,\/S) is 

proportional to J(1 - So for 

y -^tot -^tot 

large total energy much higher than production thresh¬ 
old, the cross section has mild dependence on S. Relative 
error by taking it constant is of order . 

-^tot 

Production cross section is defined as the sum of full 
scattering cross sections of all different decay channels to 
detectable stable particles. Thus, it is a detection depen¬ 
dent concept and so is our wave function renormalization 
factor, Z' which is defined as 

/•‘S'max 

Z'iE,,P,A) = J2 dSp^iS) (3) 

AeA 

where A is a set of relevent decay modes. If iftot is taken 
to be infinity, Z'(oo; all) = 1 for n > 1. 

Considering a stable particle theory helps better un¬ 
derstanding eq.(l) and (3). If we take (() to be a stable 
particle and A to be an asymptotic one particle state, 
p\ = Z5{p^ — m^). Plugging this into eq.(l), a result 
from LSZ is recovered. With the same kinds of inter¬ 
actions we are discussing in this letter, Z of the stable 
particle also vanishes as the (well defined) physical mass 
approaches one of threshold masses. On the contrary, 
from eq.(l), we see that the production of the particles 
does not change if off-shell decay modes are also con¬ 
sidered. For the stable particles, what is suppressed is 
the production of on-shell state which remains even after 
infinite time. 



It is natural to define a branching ratio to a channel A 
as 




(4) 


When we consider only the decay channels whose thresh¬ 
old masses are much below the resonance peak as de¬ 
tectable decay modes, it is easy to see eq.(4) agrees with 
the conventional definition of branching ratio. 

In total, our NWA is 


CT(ini —>■1,2,--- , n, A) ~ Z'crlim —>■ mi, - - - , m„, m\) BrA. 


For the decay channels with threshold masses much lower 
than the physical mass of the unstable particle, m\S cor¬ 
respond to the physical mass of the unstable particle. 
Scalar trilinear as an example We take scalar tri- 
linear interactions as a simple example motivated by the 
setup explaining electroweak symmetry breaking from ra¬ 
diative corrections given in [17]. The real scalar field h 
that can decay to a light real scalar b has a scalar trilinear 
interaction with another real scalar s, 


A£ = —Ah{x)s{x)"^. (6) 


where A is a dimensionful coupling and h and s have 
physical mass of mph and renormalized mass of mg, re¬ 
spectively. The physical mass is defined as a zero of 
Re(G'“^) by absorbing some parts of the self energy at 
the given scale. The remaining quantum correction to 
the inverse propagator is given by 


|2Ap 

IfiTT^ 



Wg) -Re(/(mph,mg))) 


-bimphFb, (7) 


where 

/(pAm*)=-^Axlog(AAliA±iA^) (8, 

[ 1 - ifp2 < 0 

= J 1 — |/3| tan“^ i if 0 < < 4m^ , 

[ 1-i|^|(ln^^-i7r) ifp2>4m2 

and we approximate the imaginary part from the inter¬ 
action between h and b as F^ term in eq.(7). 

In order to illustrate the influence of the threshold near 
the resonance, we focus on the limit of threshold at the 
resonance, mph = 2mg, which is shown in the second plot 
of Fig. 2. Though this example deals with the decay into 
two identical particles, most of physics discussed here is 
not affected when the decay products are two different 
particles as long as their sum is near the resonance. 

As seen in Fig. 2, the spectral density shows asym¬ 
metric cusp behavior. This distribution was discussed in 
[14] to understand the Trry system. This effect is visible 
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FIG. 2. Signal cross sections for each channels are described by spectral density p\. Black curves represent phbij?)- As the 
threshold goes to the physical mass, the shape becomes narrower and narrower. The peak position is near the physical mass 
and its width can be approximated by Fi in eq.(lO). Blue(grey) curves show pss{p^) and it is zero at and below the threshold. 
Its width is as large as r 2 in eq.(ll) and the peak position is near the threshold. The parameters are fixed by mph = 126 GeV, 
A = 0.6mph, Fi, = 4 MeV for all the three plots. These numbers are used just for an example and corresponding e in eq.(9) is 
0.155. This parameter set is also used in Fig. 3 and Fig. 4. 2ms — mph = —Fj,, 0, Fj, are used from the left. 


enough when the interaction of the resonance with the 
threshold particle is stronger than the remaining ones at 
the resonance. More specifically, the following condition 
is assumed for discussions, 


"^phr;, 

r2 

s 


< 1 , 


(9) 


where Fg = |Ap/(87rTOph) is the asymptotic partial decay 
width without phase space suppression. For e ^ 1, there 
appears two interesting parameters Fi and F 2 . 


Fi = eFf,, (10) 

F 2 = ^Ff,. (11) 

Unlike the usual interpretation of decay width coming 
from the imaginary part of the complex pole [15], these 
correspond to real deviations of two complex poles of the 
propagator in the second Riemann sheet from ruph. 

Since the spectral density is — ;^Im(G'), its behavior can 
be understood by looking at which term in the inverse 
of dressed propagator dominates as we vary energy just 
above the resonance where 


G ^ ~ mph|/3p+fTOphFg/3 + imphFt,. (12) 

: iTTiphFf, in the nearest of the resonance(mph,m.ph + Fi), 
iniphTsP near the resonance(mph + Fi, mph + F 2 ), and 
— mph away from the resonance(mph + F 2 ,oo). Thus 
the half-width of bb resonance can be represented by Fi. 
And the cross section to ss starts to appear from the 
threshold which is at the resonance and lasts up to F 2 
which corresponds to the half-width of ss resonance. 

The first (third) plot of Fig. 2 shows the spectral den¬ 
sities with the threshold mass below (above) the physi¬ 
cal mass. In the first plot, the cusp point is almost the 
threshold mass where the imaginary part is minimized. 
In the third plot, the spectral density to bb is quite sim¬ 
ilar to BW distribution except that it is narrower than 


the case with Fg = 0. For mph — 2mg > F 2 , the max¬ 
imum of bb channel resonance is mainly determined by 
the real part of the self energy and the BW distribution 
is recovered for bb channels. 

Fig. 3 describes Z'{oo-,bb), Z'(oo; 56, ss), and conven¬ 
tional Z as a function of the threshold mass. Although 
Etot does not have to be infinity, we use it for notational 
convenience to mean that the interval (S'minj 5'max) suf¬ 
ficiently covers the resonance region. As the threshold 
mass approaches the physical mass from above, the dot¬ 
ted line in Fig. 3, Z goes to zero. Thus, the conven¬ 
tional NWA shows a pathological behavior of the scalar 
h at the threshold, which is vanishing full scattering cross 
section involving bb channel. In contrast, the definition 
in eq.(3) is well defined even in this limit. For bb chan¬ 
nel, Z'(oo\bb) ~ e is small but non-zero at the physical 
mass and it gives correct approximation to the full scat¬ 
tering cross section combined with eq.(5). Alternatively, 
we can start from Z'(oo;55, ss) = 1 now with Brbb — e 
by including ss channel. Both scenarios explain that the 
full cross section to bb channel is suppressed due to the 
threshold near the resonance. 

It is interesting to ask which one of widths is the decay 
rate among Ff,, Fg, Fi, and F 2 . A survival probability is 
actually a concept depending on how unstable particles 
are prepared, or mass filtered [12]. Here, we consider 
an ideal situation in which produced unstable particles 
sufficiently cover the resonance region so that we can use 

Back to the limit 2mg = mph, major decay pattern can 
be understood by saperating the energy region (mph+Fi, 
Wph + F 2 ) and (mph + F 2 , 2mph) where the spectral den¬ 
sity 2Ep{E‘^) can be approximated by ^ 


^2 '' yjE-m-ph 


and (iv^ 




mph 


= )^, respectively. Roughly speak¬ 


ing, their effects on the survival probability P(t) are such 
that two characteristic patterns of non-exponential de- 
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FIG. 3. Z'(oo\bb) (black), Z'{oo\bb,ss) (blue/grey), and con¬ 
ventional Z (Dotted) as a function of the threshold mass. The 
conventional Z is ill-defined as its left limit and the right limit 
do not match at 2ms = m-ph while Z's are well-defined in all 
range. Z'{oo; bb) asymptotically reaches to Z as 2ms goes far 
above the resonance. 
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FIG. 4. Survival probability, P(t) obtained by a numerical 
way (Solid line), analytical approximation (dashed) and ex¬ 
ponential decay for comparison (dotted). 


cay appear, l/(r 2 t) at time much later than l/r 2 , and 
g-vT^t time before I/r 2 . 

Numerically exact result and analytic approximation 
of P(t) are shown in Fig. 4. The approximation is done 
by summing up those two functions approximated in the 
above paragraph with exponentially suppressing envelop. 
The plot shows that the decay is non-exponential. To 
help make a contrast with exponential decay, staright 
dotted lines are introduced in the plot. The reason why 
the plot looks straight at large t is that exponential and 
power functions are difficult to distinguish for short in¬ 
terval at large t in log plot. 

Classification In this section, we list possible types of 
interactions which can cause conventional Z ill-defined 
such that the discussion given in this Letter is rele¬ 
vant. The key point is from eq.( 8 ) that the leading or¬ 
der expansion of self energy around ~ i® 

/3 = Thus if mph is accidentally around 

rrii, then one loop correction overcomes the tree level 
part, — rUpj^ ~ . The leading power of /3 in 

the expansion of one loop correction determines the be¬ 
havior. If the leading power is lower than one, then Z 
vanishes as physical mass and threshold get closer and 
the ill-defined Z should be cured with the method sug¬ 
gested in this Letter. 

First, in scalar trilinear interaction, even if all the three 
scalar fields are different species, the leading order expan¬ 
sion is linear in /3. For Yukawa interaction, among 

the scalar S and fermions ipt, only the self energy of 
fermion field has a term linear in /3, while the self energy 
of scalar has (3^ due to extra factor from spinor. If the 
Yukawa interaction involves 75 , such that S''(/'i 75 ' 02 ) then 
only the self energy of (pseudo-)scalar field has the prob¬ 
lematic expansion in which the leading power is square 
root order. 

All those three cases significantly change resonance 


shape from BW form as we discussed in this letter. We 
leave gauge interaction cases as a future work. 
Conclusion We discussed physics of particle production 
and decay when the threshold is at or near the reso¬ 
nance. The conventional wave function renormalization 
factor or the derivative of the self energy is ill-defined, 
and usual Breit-Wigner approximation fails since the self 
energy varies rapidly right at the resonance. Different 
decay modes show different resonance shapes including 
the narrowly peaked resonance and the broad threshold 
resonance. Apparently there is no way to use the nar¬ 
row width approximation due to complicated physics at 
around the resonance in this case. 

In this Letter, we showed that we can still define nar¬ 
row width approximation given in eq.(5) if (i) the pro¬ 
duction cross section does not change rapidly within the 
the resonance, (ii) the resonance width is small enough 
compared to the produced invariant mass S range, and 
(iii) the branching ratio is defined as the integral of the 
specific channel spectral density as in eq.(4). There is 
no unique decay width defining the unstable particle 
and it should be found from the shape of the resonance 
(half-width) or from the decay pattern (inverse of half¬ 
lifetime). The width is then channel dependent. The 
same is true for the resonance peak position or the mass 
of the unstable particle. 

The effort to look for unanimous definition of mass and 
width for unstable particles would be in vain as the study 
in this Letter shows. Only stable particles can appear as 
asymptotic states. Whenever pathological problems ap¬ 
pear, the full self energy of the unstable particles should 
be used instead of BW approximation. Still the gener¬ 
alized NWA works, and the threshold at the resonance 
reduces the branching ratio to the other decay channels 
below the physical mass. The influence of new states at 
the resonance or slightly above is correctly taken account 
in the generalized NWA. 
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